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ABSTRACT 

A model for the main observational characteristics of the radio emission of pulsars with well 
organized drifting subpulses is presented. We propose that drifting subpulses result from the 
modulation of the radio emission mechanism due to long-wavelength drift waves in the mag- 
netosphere. The drift waves are generated at shorter wavelengths, and their nonlinear evolution 
favours accumulation in a specific azimuthal eigenmode with an integral number, m, of nodes 
encircling the magnetic pole. The electric field of the drift waves is along the magnetic field 
lines, and this modulates the distribution for particles and hence the radio emission mech- 
anism. The ratio of the frequency of the eigenmode to the rotation frequency of the star is 
insensitive to the magnetic field strength and the period of rotation, and is of order unity. The 
period, P3, of the drifting subpulses is attributed to the mismatch between this frequency and 
the nearest harmonic of the rotation frequency of the star. 

Key words: pulsars: general: radio emission - polarization. 



1 INTRODUCTION 

The observational characteristics of drifting subpu lses in pul 



sars were studied in detail in the 1 9 70s iManchestei 



1971 



1975 



j Tavlor et all Il97lt jRankin et alJ 1 1974 [Manchester et al. 
lRMiknT^83lTi98al and have attracted considerable interest in 
recent years I Ramachandra n et al Jl2002t Ivan Leeuwen et al. 2002, 
l2003tlGupta et all2004 . Earlier ideas on the polarization-angle ro- 
tation have not been confirmed, and switches between orthogonal 
pola rizations are found to play a central role in apparent depolariza- 
tion iRamachandran et all2002h . A detailed study of PSR 0809+74 
Ivan Leeuwen et al. 2002, 2003) has clarified the interplay between 
subpulse drifting and nulling. These and other investigations have 
led to renewed interest in the interpretation of subpulse drifting. 
The basic interpretation of subpulse drift involves a pattern that 
drifts around the magnetic equator such that an individual subpulse 
reappears in the pulse window after a characteristic rotation time, 
P3. A generic model, referred to as the carousel model, involves 
the drifting pattern being identified with plasma columns spaced 
periodically around a circle that rotates relative to the stellar sur- 
face about the magnetic pole. The number of columns and the pe- 
riod of the patte rn rotation have been determin ed for several pul- 
sars JPeshpande & Rankinlll999L EoOll Ivan Leeuwen et alJl20ol 
2003), and this number is of order 20. The patterns implicit in 
the carousel model may be interpreted in a variety of ways. One 
interpretation is that there are actual plasma columns, for exam- 
ple formed through sparkin g on the stellar su rface, as envisaged in 
the model of Ruderman & Sutherland 1 1975) and reconsidered re- 
cently by Gil, Melikidze & Geppert (2002). Another interpretation 



is in te rms of a magnetic field with a high order multipole com- 
ponent lAsseo &' Khechinashvili 2002). More generally, the distri- 
bution of plasma, magnetic field and other relevant quantities at a 
particular radial distance from the star can be expressed as a sum 
of spherical harmonic components, £, m, with £ = 0, 1, 2. . . ., and 
— £ ^ m ^ I, and an interpretation of the carousel model is that a 
specific such harmonic with m ~ 20 dominates all other and leads 
to the inferred pattern of m nodes or anti-nodes around a circle 
at fixed magnetic latitude. In one specific such model the pattern 
is interpreted in terms of resonant cavity modes 1 Young 2004) in a 
model in which the polar-cap region acts as the cavity. Another pos- 
sible class of models involves interpreting the spherical harmonic 
components as long-wavelength plasma wave modes such that the 
ratio of the wavelength is not negligible in comparison with the ra- 
dius of the star. This is familiar in helioseismology, in which the 
waves are sound-like modes within the star; the spectrum of vari- 
ations in the observed radiation from a star has sharp peaks that 
are associated with the natural modes of the star, described by spe- 
cific/, m values. Such a model was proposed bv lClemens & Rosea 
(2004), who invoked a s pherical mode with large £ ~ 70 and 
m = 0. A point made bv lClemens & Rosenl 120041) is that the fre- 
quency of the wave can play an important role in the interpretation 
of P3 . In the conventional carousel model the columns are assumed 
stationary in a corotating frame, and there is no counterpart of the 
frequency of oscillation of the wave. However, a model based on 
internal oscillations of the star does not lead to a natural way of 
modulating the radio emission, which is thought to originate well 
above the surface of the star, and it is unclear why the inferred pat- 
tern should be fixed relative to the magnetic axis rather than to the 
rotational axis. 
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In the present paper we propose a model for drifting subpulses 
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that is based on long- wavelength drift waves that form a period pat- 
tern around a circle cen tred on the magnetic pole. As shown by 
Kazbe ei et alJ 11991 a b), drift waves propagating across the field 
lines can be generated in the electron-positron plasma of the pulsar 
magnetosphere through a resonant instability driven by the oppo- 
sitely directed drift motions of electrons and positrons. Drift wave 
have their wave vector, k predominantly in the plane defined by 
the curvature drift, which is almost exactly perpendicular to the 
magnetic field lines. When generated, the waves have a relatively 
short wavelength compared with the dimensions of the magneto- 
sphere. Drift waves have a low group velocity, so that the wave 
energy tends to build up locally rather than propagating away, and 
their growth is limited by nonlinear processes which redistribute the 
wave energy to other scales. In particular, induced scattering leads 
to an energy transfer to lower frequencies which correspond to 
longer wavelengths. In close analogy with the formation of a Lang- 
muir condensate l Goldman 1983) in an unmagnetized nonrelativis- 
tic plasma, the wave energy accumulates in the long-wavelength 
region of phase space where there is no effective dissipation mech- 
anism and where the propagation speed of the wave energy is very 
slow. For wavelengths of order the transverse scale of magneto- 
sphere, the treatment of the drift waves as plane waves becomes 
inappropriate, and the geometry of the magnetosphere needs to 
be taken into account. A treatment in terms of spherical harmonic 
eigenmodes with specific i, m is then appropriate: an assumption 
of strict spherical geometry would be artificial, but the description 
of the azimuthal distribution in terms of eigenmodes varying as 
exp(im0) should remain valid in a more realistic model, and it is 
only the m values that are directly relevant in the model proposed 
here. Long-wavelength drift waves are then represented in terms of 
eigenmodes with specific azimuthal wave numbers corresponding 
to an integer number, m, of nodes around a circle centred on the 
magnetic pole. We argue that the nonlinear processes that transfer 
the waves to long wavelengths should result in a specific m value 
dominating. The frequency of these drift waves can be of the order 
of pulsar rotation frequency, Q. = 2n / P, where P is the pulsar pe- 
riod. Then, as in the model o f|ci emens & Rosed 120041) . the time 
P3 is attributed to a mismatch between the wave frequency and a 
harmonic of the rotation frequency. 

There are several attractive feature of the drift-wave model for 
drifting subpulses. One attractive feature is that drift waves pro- 
vide a natural way of modulating the radio emission. Specifically, 
the electric field in the waves is primarily along the magnetic field, 
and this modulates the acceleration of particles and the form of 
their distribution function. Although there is no consensus on the 
pulsar radio emission mechanism, the most widely favoured mech- 
anisms involve waves growing due to a Cerenkov instability driven 
by free energy associated with an increasing portion of the parti- 
cle momentum distribution function. This function, and hence the 
instability, is modulated by the electric field in the drift waves. At 
a specific height and magnetic latitude, the emission should then 
have m peaks in azimuth, and should vary temporally at the wave 
frequency, so that P3 results from the mismatch of this frequency 
and the nearest harmonic of the rotation frequency. Another attrac- 
tive feature is that drift waves are intrinsic to the magnetosphere. 
The pattern of drift waves, described by the distribution among the 
spherical modes I, m, is likely to change slowly with radial dis- 
tance from the star, with the drift pattern at a specific observational 
frequency determined by the dominant I, m mode at the radius at 
which that frequency is emitted (assuming some form of radius- 
to-frequency mapping). This avoids the constraint, implic it in the 
models based on sparking ( Rud erman & Sutherlandfl973) or stel- 



lar oscillations IClemens & Roserl2004) . that the pattern of drifting 
subpulses be tied to surface phenomena on the star, for which there 
is no direct evidence. 

The paper is organized as follows. The generation of the drift 
waves and their nonlinear evolution are discussed in section 2. Our 
proposed model for drifting subpulses is presented in section 3. The 
conclusions are summarized in section 4. 



2 GENERATION OF DRIFT WAVES 

Drift waves grow due to an instability driven by the curvature drift, 
which is proportional to the Lorentz factor of the particles. Al- 
though the secondary pairs are much more numerous than the beam 
particles, they have much smaller Lorentz factors. Both the disper- 
sion of the drift waves and their growth is determined by the species 
of particle with the highest curvature drift rate, which are the beam 
or 'primary' particles. 



2.1 Curvature drift instability 

Particles moving along the curved magnetic field undergo drift or- 
thogonal to the plane of field curvature. It is convenient to define a 
local Cartesian coordinate system, with the z axis along the direc- 
tion of the magnetic field and the curvature drift velocity along the 
T-axis. The drift speed is given by 



Ub — 



cyv z 

LO„Rc 



(1) 



where ujb = eB /mc, R c is the radius of curvature of the magnetic 
field line, 7 is the Lorentz factor of the particle, and v z ~ c is the 
particle velocity along the magnetic field line. 

Curvature drift can cause low-frequency, ne arly transverse 
drift waves to grow iKazbeei etailll991al 1996). These waves 
propagate across the magnetic field, so that the angle 6 between 
k and B is close to 7r/2, corresponding to k±/k z ^> 1, with 
k± = {k x + kl) 1/2 . Assuming r y(cu/iu B ) < 1, (u b /c) 2 < 
1, k z /k x <C 1 and k y —> 0. The dispersion equatio n for t he 
drift wave may be approximated by ( Kazb eei et alll991bl Fl996) 



id = k x u b + k z v z + iF, 



(2) 



with Ub given by and where T corresponds to the growth rate. 
The growth rate is largest for 



(3) 



where io v is the plasma frequency. The maximum growth rate is 
given approximately by 



1 {'— 1 k xUb . 



(4) 



The growing drift waves draw their energy from the motion 
of the beam particles along the magnetic field lines. The resonance 
condition that allows their growth requires k x Ub 7^ 0. That is, the 
drift motion of the beam particles is essential for the resonance con- 
dition to be satisfied, but the drift motion does not provide the free 
energy that drives the instability. Assuming l/k x is the order mag- 
nitude of the transverse dimension of the magnetosphere (cf. eq 7), 
from (1) one has k x u b ~ c 2 ^r /uj B ,Rl ~ 0.053(r/10 8 cm) s" 1 , 
where 7 ~ 10 7 , R t = 10 6 cm is the star's radius, r is the emis- 
sion radius, and lob* is the cyclotron frequency with the magnetic 
field B» taken on the polar cap. The growth rate of the drift wave 
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is rather small, i.e. generally r < fl, but even a small growth rate 
can be effective because the waves propagate nearly perpendicu- 
lar to the magnetic field and can stay in resonance for a long time 
3> 2n/Q. Each particle gives up only a small part of its energy 
to the waves as it propagates through an interaction region where 
the waves are growing. New particles continually pass through this 
interaction region where they spend a relatively short time so that 
the back reaction of the waves upon the particles is unimportant. 
Thus the energy accumulates in the waves without quasi-linear sat- 
uration. The wave amplitude grows until nonlinear processes redis- 
tribute the wave energy. 



2.2 Nonlinear evolution of drift waves 

Wave-wave interactions include three-wave and four-wave cou- 
plings and induced scattering. A three-wave interaction involves the 
quadratically nonlinear current which is an odd function of the sign 
of the charge, oc e 3 , so that electrons and positrons contribute with 
opposite sign, making this effect relatively weak in a pulsar plasma 
iLuo & Melrosd ll994). Four-wave interaction, which involve the 
cubic nonlinearity, oc e 4 , are intrinsically small, and tend to be im- 
portant only when three-wave processes are weak or forbidden. In- 
duced scattering may be regarded as a nonlinear absorption process, 
absorbing the beam between the scattered and unscattered waves. 
This leads to a transfer of waves from higher to lower frequencies, 
with the rate proportional to the product of the wave energies in the 
scattered and unscattered waves. Induced scattering should be the 
dominant nonlinear process in the present case. 

Induced scattering occurs when the resonant condition 



u) — <J — (k z — k' z )v z 



0. 



(5) 



is satisfied. Note that induced scattering is dominated by the pairs, 
and not by the primary particles, so that we may estimate \v z \ ~ 
1 — l/2jp. Substituting J2j into ^5) one obtains the ratio 



k x k x 
kz kL 



c 1 
u b 272 ' 



(6) 



We estimate this ratio be be in the range 10-100. The positive sign 
of k z corresponds to waves propagating from the pulsar to the light 
cylinder. 

This induced scattering transfers waves from higher to slightly 
lower frequencies, in a way analogous to the the transfer of 
Langmuir waves in an isotropic thermal plasma l TsYtovic r]| 19701 : 

sell 



Saed eev & Galee\il969tlGoldmaJl983tlMelrosJl986ir The drift 
waves generated by the instability have k z /k x <C 1, and the nonlin- 
ear transfer is from larger to smaller \k x \. From {6j it follows that 
during each scattering process the changes in the wave numbers 
satisfy [Afca;| 3> |Afe z |. This implies that the transfer of energy 
through the scales should be determined by k x alone and that the 
process is closely analogous to the induced scattering of Langmuir 
waves in an isotropic thermal plasma. As the transfer to smaller k x 
the waves move into a range of k x where there is no effective wave 
dissipation, due to the resonance condition not being satisfied for 
any particles. This energy transfer continues until the length scale 
approaches the transverse dimension, d, of the magnetosphere cor- 
responding to a minimum fcjf™ ~ 1/d. We identify d are the radius 
of a circle about the magnetic pole, and express it in terms of the 
emission radius, r, which we write as a fraction, £ <C 1, of the light 
cylinder radius, Rlc- This gives 



!/2 /?lc = (1.6 x 10 s cm) ( X 



3/2 



T~ ■ (7) 



For such small k x , the treatment in terms of plane waves becomes 
inappropriate. 

When the wavelength of the drift waves in the azimuthal direc- 
tion is of order d, the waves need to be treated in terms of azimuthal 
harmonics, oc exp(im</>), rather than plane waves. Qualitatively, 
the important point is that relevant waves have an integral number, 
m, of nodes and anti-nodes around the circumference of the circle. 
The wave energy accumulates in this wavelength region provided 
that the rate of nonlinear pumping into the region with k x ~ fc™ ln 
is sufficiently high. Specifically, the nonlinear transfer must occur 
at a faster rate than the drift waves propagate out of the region, al- 
lowing the 'reservoir' of energy (in wave number space) at k x » 
10 _s cm" 1 to be maintained. The propagation time for waves with 



k z /k x w 10" 



(£/n)(k/*») w 2(C/0.1)(P/1 s) s. Con- 



sequently, the reservoir is exhausted in time r. We need to com- 
pare this with the characteristic time for the nonlinear transfer to 

h ^ £- min 
l^x ^ x 

The nonlinear transfer rate due to induced scattering may be 
estimated as 



: 2> 



(8) 



where r\ = U /n p m e c 2 -yp is the ratio of the energy density of the 
waves, U, to the energy density in the beam particles. It follows 
from ^8) that setting Fnl equal to the linear growth rate implies 
only a modest energy density in the waves, r\ <C 1. It then follows 
that the transfer time, 1/Tnl, for the drift waves to small k x is of 
the order of the time r. We assume that a balance is established 
between input of wave energy due to the instability at relatively 
large k x and nonlinear transfer of energy to smaller k x . 

We conclude that it is plausible that drift waves do accumu- 
late in the long-wavelength region. The nonlinear transfer should 
favour one specific eigenmode becoming dominant, corresponding 
to a specific azimuthal wave pattern, that is, a pattern with a specific 
integer number, m, of wavelengths around a circle centered on the 
magnetic pole. 



2.3 Frequency of the drift eigenmodes 

A remarkable feature of this model for the generation of low- 
frequency drift waves is that the resulting frequency of the waves 
is of the same order of magnitude as the rotation frequency of the 
star. In particular, the ratio of the wave frequency to the rotation 
frequency depends only on the mode number, m, and the ratio, 
£ = r/Ri,c, of the emission radius to the light cylinder radius. 
This occurs because the Lorentz factor of the beam particles is lim- 
ited by the potential available, which is proportional to B»f2 2 , so 
that the curvature drift speed Q is independent of the strength of 
the magnetic field, B, . Let 7 max be the maximum Lorentz factor, 
corresponding to the maximum potential energy available (between 
the magnetic axis and the last closed field lines) in a vacuum model. 
One has 7 max ~ £l 2 UB*R'i / (? , and with R c w £ 1//2 i?LC this im- 
plies 



If* ~ _ 

C ""/max 



(9) 



Then for fc™ m = m/d with d given by 0, the ratio of the frequen- 
cies becomes 



IT 



7 



(10) 
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The factor 7 /7 max for the primary particles depends on the details 
of the model. Nevertheless, it is evident that for plausible values of 
m ~ 20 and £ ~ 0.1 the two frequencies are of the same order 
of magnitude. This is an important feature of our model in that it 
provides a natural explanation of the period, P3, associated with 
drifting subpulses, as discussed in the next section. 

The estimate 1 1 01 is derived assuming that the term k z v z ~ 
k z c is no larger than fc™ m itf,. A (relativistic) particle is accelerated 
along the magnetic field due to the electric field of the drift wave, 
5E\\, and k z c < uj m implies that this acceleration is temporally 
limited rather than spatially limited. 



2.4 Amplitude of the drift waves 

The effect of the drift waves on the particle distribution that is re- 
sponsible for the radio emission may be estimated by considering 
the amplitude, £7, of the oscillating part of the Lorentz factor of an 
electron or a positron due to the forced motion in the wave. Let SB 
be the magnetic amplitude of the wave. Then the relation between 
^7 and SB may be determined as follows. 

The oscillations of a particle in a wave may be either approx- 
imately spatial or approximately temporal, depending on whether 
the phase velocity of the wave is much less or much greater, re- 
spectively, than the velocity of the particle. In the present case, the 
relevant motion is along the magnetic field lines and the relevant 
phase velocity is ui m /k z . The dispersion relation u m = k m Ub and 
the assumption k z c <C lo then imply 



8E« = lib SB, 



S'y ' 



e5E« 



uj m m e c 



(11) 



where u) v is the plasma frequency for the pairs. The amplitude SB 
may be expressed in terms of the ratio r\, introduced in {SJ. One has 



n - 



SB 2 



SB 2 oj% 



/i n p m e c^7p 



B 2 



(12) 



where lo v is the plasma frequency for the pairs. The amplitude of 
the oscillations in the Lorentz factor may then be written in the 
form 



^7 



1/2 



LUp Ub 



(13) 



with M(,/c given by The assumption that the waves have small 
amplitudes requires not only ?7 <JC 1 but also S^y/fp < 1. The drift 
waves modulate the radio emission only if Sy/jp is large enough 
to modify the emission process, as discussed below. 



3 DRIFTING SUBPULSE PHENOMENON 

A model for subpulse drift based on long-wavelength drift waves 
has a number of attractive features. As already noted, the waves 
should accumulate to form a specific azimuthal pattern, qualita- 
tively consistent with the carousel model. Moreover, the fact that 
the frequency of the waves is comparable with the rotation fre- 
quency of the star implies that this pattern drifts in azimuth, again 
consistent with the observed drifting patterns, and providing a natu- 
ral explanation for P3. Another favourable feature of the drift- wave 
model is that the electric and magnetic fields of the waves provide 
a natural mechanism for modulating the radio emission. 



3.1 Modulation of the emission process 

Any model for pulsar radio emission must involve so me form 
of co herent emission, and most emission mechanisms 1 Melrose 
1 19951) involve a maser-type or reactive instability. The most widely 
favoured emission processes involve a beam instability generat- 
ing waves that indirectly produce the escaping radiation. There 
are al so several maser emis s ion pr ocesses: maser curvature emis- 
sion JLuo & Melrose! ll992 Ml995l) . linear acceleration emission 
( Melroje]jlj7jyRowej 1199511 and free-electron maser emission 
(Funs & Kuiioers 2004). A common feature of all such instabil- 
ity mechanisms is that they are driven by free energy associated 
with a peak in the momentum distribution. Specifically, if f(p) is 
the one-dimensional momentum distribution function of the par- 
ticles, these instabilities are driven by particles with df(p)/dp > 
0. Such a rising portion of the distribution function is confined 
to the lowest energies, below the maximum in f(p). Typically 
the growth rate for these instabilities decreases with increasing 
Lorentz factor at the peak in the distribution function. Recent es- 
timates suggest a relatively modest value for the peak Lorentz fac- 
tor, 7p ~ 10-100 Izhang & HardiiI3l200ri li-iibschman & Aronsl 
l200lilArendt & EilekfcOOa) . A modulation Sy/y p ~ 0.1-1 would 
strongly modify the growth of the waves in any such instability. 

The magnetic field in the drift waves also affects the pattern of 
any radio emission. The magnetic field in the waves is nearly per- 
pendicular to the pulsar magnetic field B, and leads to as a periodic 
perturbation of it. This perturbation affects the angle of emission, 
which is centred on the direction of the magnetic field, causing the 
emission angle relative to the star to oscillate in the polar direction 
away from the magnetic axis. Such oscillations imply that the emis- 
sion pattern from dipolar field lines originating from a fixed mag- 
netic latitude on the star form a pattern with m oscillations in the 
polar direction about a circle on the sky around the magnetic pole. 
The angle through which the direction of emission is deflected is 
~ SB IB, cf. (T2T 

In principle, the oscillations in the magnetic field, SB, 
might also act as a wig gler field in a free-electron maser model 
iFung & Kuiioers 2004). However, the frequency of the emission 
would be ~ k x cy 2 , which is too low to be of observational interest 
from k x ~ fc™ m and 7 P ~ 10-100 for the bulk of the pairs. 

3.2 Interpretation of P3 

As already noted, the frequency, u m , of the drift waves is of order 
the rotation frequency, Q, and this provides an interpretation of the 
period P3 in terms of a beat between these frequencies. Suppose 
that u) m is close to an integer, n, time the rotation frequency, 



ifl ± A. 



(14) 



Then the drifting patter n repeats with a period P3 = 2tt/A 
IClemens & RoserJliool) . For A > the subpulse drift is in the 
direction of the pulsar rotation and for A < 0, the wave phase 
falls behind the pulsar rotation and subpulse drift is in the oppo- 
site direction. After the time P3 = 2ir/|A| the phase returns to the 
original position and the pattern repeats. 

The ratio of the numbers n and m may be estimated as fol- 
lows. The wave number is fc™ m = m/d, and the frequency is 
uj m = mub /d. Assuming a dipolar field, this give 
n 7 

m 7ma: 

where j5J and J I4l are used. The ratio i 151 for emission at a height 
corresponding to £ ~ 0.1 and for 7/7 max ~ 1 is consistent with 



(15) 
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n of a few and m ~ 20. Such numbers can plausibly explain the 
main features of subpulse drifting. 



3.3 Further implications of the model 

Our model differs from all other models for subpulse drifting in 
that it relies purely on magnetospheric oscillations and not on ef- 
fects related to the surface of the star. In our discussion we have 
assumed that the emission height is a fraction £ ~ 0.1 of the dis- 
tance to the light cylinder. Although the results are not particularly 
sensitive to this choice of £, it is clear that the model cannot ex- 
plain radio emission originating from close to the surface of the 
star, e.g. Kiiak & Gil 1 1998). The purely magnetospheric charac- 
ter of the modulation of the radio emission is an attractive feature 
of the model, but it does imply that the source of the radio emis- 
sion is relatively high in the magnetosphere. In particular, it is not 
consistent with some rec ent models for radio emission from near 
the surface of the star, cf. lGil. Lvubarskv & Melikidze 1 2004J). The 
model also provides scope for interpreting a variety of phenomena. 
We comment briefly on some possibilities. 

Pulsars PSR 2303+30 and PSR 1944+17 exhibit a complicated 
picture of drifting subpulses that is usually explained by the exis- 
tence of two or more drift systems l Manchester & Tavlor 1977). 
An interpretation within the framework of the present model is 
that two different azimuthal harmonics (two different m values) 
for the drift waves are excited with approximately the same ampli- 
tude. The phenomenon of nulling is associated with drifting sub- 
pulses, and one possible interpretation is that emission is possible 
only in the presence of the drift waves (occurring when the oscil- 
lations in 7 P have their minimum value), and that for some rea- 
son the pattern of drift waves disappears. The transient features 
observed as the drifting subpulses are re-established following a 
null l van Leeuwen et al. 2003) may be interpreted in terms of the 
re-establishment of the drift-wave pattern. Although these sugges- 
tions are only speculative, the important qualitative point is that the 
drift-wave model provides a purely magnetospheric framework for 
interpreting phenomena associated with drifting subpulses, without 
any reference to surface phenomena on the star. 



4 CONCLUSIONS 

In this paper we attribute the subpulse drifting phenomenon to mod- 
ulation of the radio wave generation region by drift waves. Rel- 
atively small scale drift waves are generated by a beam instabil- 
ity due to the curvature drift of 'primary' particles. Nonlinear pro- 
cesses, specifically induced scattering of the drift waves by plasma 
particles, transfers the drift wave energy to longer wavelengths. At 
sufficiently long wavelengths the geometry of the magnetosphere 
becomes important, and the eigenmodes correspond to specific az- 
imuthal numbers m. We argue that the nonlinear wave processes 
should favour one eigenmode becoming dominant, so that the wave 
energy accumulate in waves with this m. The waves then form a 
pattern with an integral number, m, of nodes around a circle cen- 
tre on the magnetic pole. We suggest that the m nodes or anti- 
nodes around a circle correspond to the number of subsources in the 
carousel model. This pattern varies in time at the wave frequency, 
and the mismatch, A = uo — nQ, between the wave frequency 
and the nearest harmonic of the rotation frequency determines the 
period, P3 = 2ir/ A, at which the pattern repeats. For the model 
to account for the main features of the observed subpulse drift re- 



quires that n = lj/Q be of order unity, and m be of order 20. For 
plausible parameters we find m/n of the required order, cf. 1151 . 

These long-wavelength drift waves can affect the radio wave 
generation process by modulating the low-energy portion of the 
particle distribution function due to acceleration along the magnetic 
field lines by the electric field of the drift waves. Most favoured 
emission processes involve an instability driven by the rising por- 
tion, df(p)/dp > 0, of the distribution function, and the emission 
processes are sensitive to modifications in this low-energy portion 
of the distribution function. The magnetic field in the drift waves 
also modifies the local direction of the magnetic field lines, leading 
to a change in the angle of emission relative to the star. 

Qualitatively, perhaps the most important feature of the drift 
wave model is that it is purely magnetospheric, and is not related to 
processes on the surface of the star. We envisage emission from a 
height that is a fraction, £, of the light-cylinder radius, with £ ~ 0. 1 
in our numerical estimates, and the drift-wave pattern is assumed 
to be set up at this height. Our favored suggestion for modulation 
of the emission process, to produce drifting subpulses, is that the 
parallel electric field in the drift waves modifies the distribution 
of particles, and allows maser emission to be effective only in m 
localized regions, implying a carousel-like model. Another possi- 
bility is that the effect of the drift waves is to modulate the angle 
of emission so that it oscillates in and out of sight as the carousel 
rotates. 
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